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I. INTRODUCTION 

The Farey fraction spin chain, which we study here, is one example 
of a set of closely related one-dimensional models (see [l| and Q for de- 
tails) which are of interest in both statistical mechanics and number theory. 
From the statistical mechanics point of view, there are a number of results, 
mainly for h — 0. All the models have the same free energy f((3,h = 0), 
and exhibit a phase transition, at finite temperature, that is rather unusual. 
It is known, rigorously, to lie on the border between first- and second-order 
with the asymptotic form /(/?, 0) ~ t/logt (where t = 1 — (3/P c )- In the 
low-temperature state /(/3, 0) = and the magnetization is saturated, so 
there are no thermal effects at all. For (3 < (3 C , f(f3, 0) < and the magne- 
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tization vanishes (see [3| for details). In [jj it is proven that the saturated 
state persists for h ^ and (3 > (3 C ; when h > 0, the magnetization m = 1 
and when h < 0, m = — 1. The interesting question is how these states 
relate to the high-temperature state. In this question is addressed via 
a renormalization group calculation which finds, among other results, a 
phase diagram that is the same as illustrated in Figure below. However, 
these models have long-range interactions, so the applicability of the renor- 
malization group might be questioned. It was in fact the desire to verify 
the results at non-zero h in [jj that motivated this work (the results of 1] 
for h = agree with the rigorous behavior) . 

The Farey fraction spin chain has also led to some new results in number 
theory 0, IE 01 ■ Additionally, as explained below, there is a connection to 
dynamical systems-in fact, the asymptotic form of /(/?, 0) mentioned comes 
from a result for dynamical systems Furthermore, chaotic behavior is 
even exhibited by certain statistical quantities. In particular, [f| proves 
that the "density of states" for the infinite chain does not exist-it is a 
distribution. 

In what follows, we analyze the Farey fraction spin chain in a particular 
approximation, in which the energy of a configuration is described by single 
cluster energies. Within this approximation, the model becomes exactly 
solvable (and is closely related to the "necklace" or "bead" models of Fisher 
and Felderhof |8||). Furthermore, our results for h ^ agree, in the main, 
with a previous analysis Q that makes use of mean field theory and the 
renormalization group. There are, however, some intriguing differences. 

The Farey fraction spin chain |3, |9( may be defined as a chain of N spins 
ai; i — 1,2, ... N with two possible states Ui € {|, J,}. Using the matrices 

M = (\ f) and A l = (o f) ' (1) 
we define the energy of a configuration of N spins as 

N N 

E N ({a i },h) = log(Trl[A !7i )-hJ2(XT(<ri)-Xl(<7i)), (2) 

where XT( CT i) = 1 (0) f° r a i = T (I) so that it counts the number of up spins; 
Xli&i) i s defined similarly to count the number of down spins. The cyclic 
invariance of the trace in makes the system translationally invariant. 
The partition function is then given as a sum over 2 N spin configuration 

Z N {P,h) = £V/MM{*«},fc) . ( 3 ) 

Our focus is the limiting free energy 

-0mh)=Um ±-log Z N (/3, h) . (4) 
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The paper is organized as follows. In section^ we utilize the thermody- 
namic formalism to give a dynamical systems interpretation of the Farey 
fraction spin chain and describe how this connection may be used to ob- 
tain an effective cluster approximation. The analysis of the Farey fraction 
spin chain within the cluster approximation is then described in section lTnl 
leading to explicit equations for the free energy and the phase boundaries. 
Section lTVI contains the calculation of the scaling properties near the critical 
point. A summary and comparison with the results from Q are contained 
in section[V] Appendix lAl contains the derivation of the asymptotics of the 
cluster partition function. 

The remainder of this section deals with a reformulation of the model in 
terms of clusters of consecutive spins of equal state and closes with a brief 
discussion of our strategy for calculating /(/3, h) in the cluster approxima- 
tion. 

Iteration of the matrices Af and A± leads to 

AS = (I f] and A? = ( I ?V (5) 



n lj i \0 l / 

One notices that while some matrix elements increase in size, the zero 
field energy for the associated configurations remains constant, as TrA™ = 
TrA™ = 2. These two states, in fact, are the ground states at zero field. 
(They are also responsible for the low temperature thermodynamics @.) 

The energy is increased considerably, however, once a change of spin oc- 
curs. It therefore is useful to think of a general configuration as a sequence 
of clusters of consecutive spins of equal state (irrespective of whether the 
state is j or |, as in zero field the energy is invariant under spin flip). If 
there is no change of spin at all, one has (as mentioned) 2V(]X_ x A ai ) = 2. 
Once there is a change of spin, we can take advantage of the cylic invariance 
of the trace to make the configuration begin with Aj (resp. A±) and end 
with Ai (resp. Aj). Thus the total number of spin changes 2K must be 
even (with K > 1), and we can describe such a spin configuration {(Ji}f =1 
by a sequence of 2K clusters of size > 1 with X)fe=i n k — A. Therefore, 
using 

A i = SA^S- 1 with S = U J) = S- 1 (6) 
we can write for any configuration with K > 1 

N 2K /n 1 \ 

Tr J] K t - Tr [] M rik with M n „ = A" fc S = K M ■ (7) 

i=l fe=l ^ ' 



Let us now suppose that we could find a meaningful approximation of 
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the form 

2K 2K 

Tr [] M nk « J] e e "* . (8) 

k=l k=l 

In this case, the energy simplifies considerably and (for all but the ground 
states) we can write 

2K 

ENdai}, ft) « £ ( £ « fc - {-if-'hnu) . (9) 
fc=i 

Working in the grand canonical ensemble, we show in section ITTT1 that one 
obtains from this approximation an exact expression for the limiting free 
energy f((3,h), given by 

Aie^-V , 0) A(e /3(/+ ' l) ,(3) = 1, (10) 
with the cluster generating function 

oc 

A(z,/3) = ^^e-' 3£ " . (11) 

n=l 

It follows from the approximation discussed in section ITTTI that this gener- 
ating function has a radius of convergence of 1. Hence for (3 > f3 c (h) the 
limiting free energy is f(f3,h) = —\h\, which for > /3 C (0) agrees with the 
rigorously known result in [l|. Thus the phase boundary (3 c (h) is given by 

A(l,/3 c )A( e - 2 ^l' l l,/3 c ) = l. (12) 

In particular, for zero field h = 0, the free energy /(/3, 0) and critical 
temperature /3 C (0) follow from 

A(e /5/ ,/3) = l and A(1,/3 C ) = 1, (13) 

respectively. 

II. THERMODYNAMIC FORMALISM AND CLUSTER 
APPROXIMATION 

In order to proceed further, we consider the thermodynamics of the Farey 
tree 101. This can be recast in a transfer operator formulation associated 
with the iteration of an interval map (see and By modifying this 

interval map we arrive at the desired cluster approximation. Note that the 
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Farey tree is known, rigorously, to have the same free energy /(/3, h = 0) 
as the Farey fraction spin chain (with no external field) Q . 

The Farey tree is generated by the Farey map defined on the unit interval 
[0,1], which is defined as 

\ = (1 - , if 1/2 < x < 1, v ' 

(see Figure QJ. We denote the inverses by F (x) = fo~ (x) — x/(l + x) and 







| fo n "'fi(x) 




fi(x)\ \ \\\ 

















f(x) 



FIG. 1: Farey map and first-return map on the interval [1/2, 1]. The first-return 
map is given explicitly by the branches /o" _1 /i for n G N. Their extension to all 
of [1/2, 1] is shown by dashed lines. 

F\(x) — fi~ 1 (x) — 1/(1 + x). The associated transfer operator is formally 
given by 

Cpftx) = \F , (x)\^(Fo(x)) + \F 1 \x)fcj ) (F 1 (x)) 

(15) 



{1 + xfP 



l + x 



1 + x 



Therefore, the JV-fold iterated operator Cp 4>(x) consists of 2 N terms of the 
form 

|(F T1 o F T2 o . . . o F rjv )'(a;)|' 3 0(F Tl o F T2 o . . . o F TN (x)) (16) 
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with Ti £ {0, 1}. As we are dealing with iterations of Mobius transforma- 
tions of the form ff^j with determinant ±1, we can alternatively consider 
multiplication of the associated matrices. (In a slight abuse of notation, we 
shall denote the Mobius transformation and the associated matrix by the 
same symbol.) We find for instance 

/#i(o)=x;<*S}> (it) 

{Ti} 

where c?{ n } is just the bottom right entry of the matrix product 

J[F Ti where F = Q and F x = ^ ^ . (18) 

Now 

iff- 1 *! = (? i) = m « ( ig ) 

which immediately suggests a cluster approximation analogous to that for 
the Farey model discussed in section [I] 

There are differences with the Farey fraction spin chain, however. One 
of them concerns details in the allowed clusterings. Both the Farey fraction 
spin chain partition function and the ./V-fold iterated transfer operator are 
expressible in sums containing 2 N terms. However, in the clustering repre- 
sentation of the Farey fraction spin chain one needs to take the cyclicity of 
the trace into account. This means that in the Farey fraction spin chain all 
excited states contain an even number of clusters. These issues (including 
multiplicity) are dealt with in section ITTT1 

More importantly, the energies are different. Note that the matrices 
Fq and Fi used to generate the d{ Ti } i n 03 differ from the matrices 
and A± employed in the Farey fraction spin chain. Therefore d{ T .y in 117|) 
cannot be compared directly with the trace (which is given in a similar 
way by a^ ai j + in Despite this, both expressions lead to the 

same free energy f([3,h — 0). This was already known via the argument 
in that the largest eigenvalue of Cp is e 2/3 ^ 2/3 \ and the corresponding 
eigenvector is positive. More recently, the connection has been shown to be 
more direct. || proves that (see fT^) 1(0) = 2Z%_ 1 (2/3), where Z§ is 
the partition function of the Knauf spin chain, which is rigorously known 
to have the same free energy f(fi,h= 0) as the Farey fraction spin chain 
(with no external field). Therefore it is reasonable to use the Farey tree 
cluster energies (see (1241) ) in a cluster approximation for the Farey fraction 
spin chain. 

We now come to the main point of this section, that one can construct a 
piecewise linear version of the Farey map which captures its essential fea- 
tures while being significantly easier to analyse. This is done by linearizing 
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FIG. 2: Linearized Farey map and first-return map on the interval [1/2, 1]. The 
first-return map is given explicitly by the branches /o n_1 /i for n G N, which are 
linear maps onto [1/2, 1]. Their extension to all of [1/2, 1] is shown in dashed 
lines. 



the map between the inverse images F fc (l/2) (see Figure 0. This sequence 
tends to zero, so that the structure of the map near the fixed point at zero 
is preserved under linearization. Due to this fact the critical dynamical 
properties of the linearized map and the Farey map are still closely related 
[l3| . In fact, the spectral radius of both associated transfer operators shows 
the same type of singular behavior at /3 C [l4|. 

Next, we consider the first-return map on the interval J = [1/2, 1], which 
is the map obtained by repeated iteration of / until f n (x) lies again in J. 
This map is the key to understanding intermittency in the Farey map |15| . 
Furthermore, it is responsible for making the particular piecewise linear 
map that we use easier to analyse than the original smooth one. If we define 
n(x) = min{n > l\f n (x) € J} then the first-return map can be written as 
g{x) — f n ( x '(x). Due to our particular choice of J the first-return map g 
becomes particularly simple. Its branches are given by /o™ 1 fi for nel 
(see Figure^. One can define a (suitably modified) transfer operator M z ,p 
for the first-return map, given by 

oo 

M x ,p4>ix) = z»\(F 1 F Q n - 1 )'(x)\l 3 <KF 1 F Q n - 1 (x)) . (20) 

71=1 
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As the first-return map is only denned on J, this operator acts on functions 
with domain J. However, using the conjugacy C(f>(x) = \f\ (x)\@ <j)(fi{x)) 
we obtain an equivalent conjugate operator C _1 M. Z $C acting on functions 
with domain [0, 1] and given by 

oo 

C-'M^C^x) = z n \(F Q n ~ 1 F 1 )'(x)f<j)(F Q n ~ 1 F 1 (x)) , (21) 

n=l 

which allows us to make the identification with M n = Fo n ~ 1 Fi. 

The crucial observation is that the operator spectra of Cg and M z ,/3 are 
related 0,0], m the sense that A = z^ 1 is an eigenvalue of Lg if and only 
if 1 is an eigenvalue of M. z ,p (f° r a rigorous formulation see HqV 

The important consequence of the particular linearization chosen is that 
the piecewise linearised map replaces the first-return map on J by a first- 
return map with branches Z^ 1-1 fi{x) which are linear and onto, with slopes 
n(n + 1). It follows that the eigenfunction associated with the leading 
eigenvalue A(z,(3) of M. z ,p becomes constant, and that this eigenvalue is 
given explicitly by 

oo n 

A(z,/3) = V Z . (22) 

' (n(n + 1 )) p 

Note that A(z,/J) extends to an analytic function in the complex z-plane 
cut from 1 to infinity for all (3. 

Reformulated in terms of approximations to iterates of the transfer op- 
erator Cp l|17|l. the linearization implies the approximation 

and we recognize that the energy has become a sum of cluster energies 

e„ = logn(n + l) . (24) 

In other words, the leading eigenvalue 11' I'D of the modified transfer operator 
of the linearized Farey map is identical with the cluster generating function 
obtained from (fTTJl and J21J. 

We have mapped the thermodynamics of the Farey tree to a dynamical 
system given by an interval map, and introduced a linearization of this 
interval map. In this way we obtain a cluster energy approximation to the 
original energy expression. This replaces the original model, which is very 
difficult to work with, with an approximation that is quite tractable. In 
the next section, we apply this cluster approximation to the Farey fraction 
spin chain in a field and determine its consequences. 
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III. THERMODYNAMICS IN THE CLUSTER 
APPROXIMATION 



As discussed in the introduction, the Farey partition function can be 
written in terms of the cluster representation {rifc}?=i- There are two 
ground states with all N spins either up or down. The excited states 
consist of configurations with 2K alternating clusters of spins with length 
n k > 1- This leads to a degeneracy n%, since (as discussed) we insist that 
the first and last clusters must have opposite spin directions. 

Now the first cluster of size n\ has either all spins up or all spins down. 
It is therefore convenient to split the partition function into two terms, 
according to the state of the first spin. If wc denote the partition function 
of configurations with <j\ =\ by Z^ N (P,h) and the partition function of 
configurations with =! by Z^(/3,h), then Zj^(f3,h) = zjj(/3,—h) and 
we obtain 



«{1 L«/2J} 



(25) 



The partition function can then be obtained as 

Z N {fl, h) = Z]f{p, h) + Z l N (p, h) = Zlifi, h) + Zl{p, -h) . (26) 
Using the approximation discussed in Section ITT1 leads to 



2 A' 
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(P,h)* (2e hN ) + £ nin^ Ws+( " iritet) • (27) 



k=l 



where now 



e„ = -logn(n + l), (28) 

with the factor 1/2 arising from (|17l) . (Note that the energy l|28|l of a 
cluster is, for large n, the same as the energy of the lowest excited state of 
a chain of length n.) From this, one can determine the limiting free energy 
f(P,h). Note also that this approximation does not change f3 c . This is 
because the change in spectral radius of the transfer operator occurs at the 
same temperature after linearization 0, . 

Passing to the grand canonical ensemble, we write 

oo 

G T (z,/?» = J2z N zl(J3,h). (29) 
JV=1 
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G^(z, (}, h) and G(z, /3, /i) are then defined similarly. One finds 

oo oo 2K 

Jf=l m,...,n 2 K=l fe=l 
— flfr oo cxd 

K=lm=l 

oo 

... £ {ze ph ) n2K e'^K . (30) 
As in the introduction, we define the cluster generating function 

oo 

A(z,(3) = J2z n e~ 0en ■ (31) 

n=l 

After some transformations we arrive at (note the close resemblance to the 
results in section 6 of 0-one difference being that the parameter if) is a 
function of f} here) 

rU o ^ ~ o-g ze ~ 0H ze-^d 1 A(ze-^'\f3)A(ze0\P) 
U(z,p,n)~z 1 _ ze - l3h + i-A(ze-Ph,p)A(zeP h ,f3) ' 1 ' 

This is perhaps most easily seen by expanding equation backwards. 
The denominator in the second term corresponds to the sum over n\ and 
ri2- The partial derivative leads to the multiplying factor tlx- Expanding 
the denominator into a geometric series produces the product over the sums 
with summation index rij with i > 3. We now get 

G(z,(3,h) = GHz,/3,h) + G l (z,p,h), (33) 

where G^(z,(3,h) = G' (z, /3,—h). The limiting free energy is then given as 



/3f(P,h) = log z c (0,h) , (34) 

where z c {f3, h) is the smallest singularity of G(z, (3, h) for z on the positive 
real axis. This singularity is reached at the smallest positive solution z c of 
one of the three equations 

ze- 0h = l, ze (3h = l, A(ze- f}h ,[3)A(ze (3h ,f3) = l. (35) 

The first two equations correspond to the two fully magnetized phases, 
where z c — e~P\ h \ or, in terms of the free energy, / = —\h\. The third 
equation corresponds to the high-temperature phase. In terms of the free 
energy, 



A(e^- h \p)A(e^ +h \p) = l 



(36) 
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Since the cluster generating function A(z, (3) has radius of convergence z = 
1, the phase boundary is given by / = —\h\ = — ft c , and h c (/3) is determined 

by 

A(l : (3)A(e- 2/3h %(3) = 1 . (37) 

The three phases meet at a critical point given by ft = and A(l, /3 C ) = 1, 
i.e. C = 2. 



IV. SCALING OF THE FREE ENERGY 



In this section we calculate the asymptotic behavior of the solutions 
/(/3, ft) and h c (0) of (|36|l and (|37[1 . respectively, near the critical point, 
which is given by ft — and /3 = @ c . For this we need the asymptotic be- 
havior of the cluster partition function for (3 — > f3 c and z — > 1. Introducing 
the reduced temperature t = 1 — (3/(3 c , we find 

A(z,/3)~ 1 + Cf + (1 - z) log(l - z) (38) 

(a derivation is given in Appendix IaT). Inserting z = e /3< ^ ± ' 1 ) into (|38|l gives 
1 — z ~ — /3 C (/ ± and we obtain to leading order 

A( e W ±/l \ „ i + ct - f3 c (f ± ft) log[-/3 c (/ ± ft)] . (39) 

Thus, (|36|l implies that to leading order 

2Ct ~ /3 C (/ + ft) log[-/3 c (/ + ft)] + &(/ - ft) Iog[-/3 c (-/ + ft)] . (40) 

In the field free case (ft = 0) this simplifies to 

Ct~/3 c /log(-/3 c /), (41) 

which can be inverted to give 

f ~iw t - (42) 

For the phase boundary (— / = |ft| = ft c ) we find, using (|37|) 

2Ct~ ~2/3 c ft c log(2/3 c ft c ) , (43) 
which can be inverted to give 

P c log* 

(see Figure OJl . Note that this result implies that on the phase boundary 
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FIG. 3: Phase boundaries for the Farey fraction spin chain near the critical 
point, showing the disordered high-temperature phase (top) and the two fully 
magnetized phases (bottom left and right). 



to leading order /(/3, h c {(3)) ~ f((3, h = 0) as f3 — > (3 C , and differences only 
appear in higher order terms. 

Finally, we compute the leading correction to the zero-field free energy 
in h. Assuming \h\ <C — / (this assumption of course breaks down near the 
phase boundary) in gives 

2Ct ~ 2/3 c f log(-/3 c /) + c h 2 /f . (45) 

Inverting this finally leads to 

(46 ) 

3 f3 c logt 2C t y ' 

for \h\ < |</logt|. 

Equations (|42|l . 144|l . and 14tjf> are our main results. We discuss them in 
the next section. 



V. SUMMARY AND DISCUSSION 



In this paper we have calculated the free energy f((3, h) of the Farey 
fraction spin chain in an external field h by making use of a cluster approx- 
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imation for the energy of the excited state spin configurations. 

We conclude with a summary of our results and comparison to previous 
work. 

• In the case of zero field, we find 

f(P,h = 0)~^-±- where t = 1 - 0/0 e . 

Pc log t 

The temperature dependence of this result agrees with the known 
rigorous result the renormalization group calculations Q, and a 
rigorous analysis of the (non-linearized) Farey transfer operator |ll| . 

• For the phase boundary, where h c — \h\ = — /, we find 

M/3) ~ — Jr^r where t = 1 - 0/0 c . 

Pc log t 

The temperature dependence again agrees with the renormaliza- 
tion group calculations Q. Additionally, we find that the con- 
stants in f({3, 0) and h c ((3) are equal, i.e. we have to leading order 
h c (P) ~ ~f(0,h = 0), (where h c {(3) = -f((3,h c (0))). This is be- 
yond what renormalization group calculations can predict, since there 
are several undetermined constants in that case (see pj for details). 
What is more interesting is that this equality is not consistent with 
the renormalization group results, as explained below. 

• The change of the free energy for small fields is given by 

fifi, h )-lTT^-^% where 1 = 1 " PIP- and W « l*/log*|- 
p c log t 2C t 

This is not in accordance with renormalization group calculations. In 
that case one has pj 

logt t 

i.e. a correction term of the order of h 2 logt/t, with a and b unde- 
termined constants. Setting this expression equal to — h c in order 
to determine the phase boundary, one finds that if the constants in 
/(/?, 0) and h c are equal, i.e. if the cluster results just mentioned 
hold, one must have 6 = 0. This is, in a sense, consistent, since it 
might imply that the leading correction to the free energy for finite 
h is of higher order than h 2 \ogt/t, and therefore could indeed be 
h 2 /t. However, it does not seem to be possible to alter the renor- 
malization calculation to obtain this while keeping the correct form 
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for f(/3, 0), the free energy at h = 0. Setting 6 = implies that the 
parameter x — in 1]. If one then includes a higher order term in 
the flow equation for u (equation (10) in lj), the result for /(/?, 0) 
is no longer correct. Another possibility is that u is a "dangerous 
irrelevant variable" (or, more precisely, a "dangerous marginal vari- 
able" ) , and the finite-field correction term to the free energy takes the 
form h(£o) 2 / (t(£o)\nu(£o)) far from the critical point, which results 
in a leading contribution to / in agreement with the cluster result 
found herein. Since very little is understood about u, however, such 
an assumption is completely ad hoc. So the question of the correct 
leading term for the free energy at finite h remains open. Now, one 
might question the applicability of the renormalization group to this 
model, due to the presence of long-range forces. The cluster results, 
being more closely tailored to the Farey fraction spin chain, are per- 
haps more compelling, but they are not rigorous either. Therefore, it 
would be interesting to know what the correct behavior is. A rigor- 
ous asymptotic analysis of the transfer operator for finite h appears 
possible [16J, and should answer this question. 




• As mentioned in the Introduction, there is a set of closely related 
models, including the Farey fraction spin chain, which all have the 
same free energy f(/3, h = 0) at zero external field. In addition, 
the magnetization is the same (see 3] for details). One therefore 
expects, according to scaling theory, that they have the same free 
energy for h ^ as well. The analysis in p| makes this assumption. 
However, it has not been proven. Given the presence of long-range 
forces in these systems, one might doubt its validity. The results 
obtained here do support it, in that the renormalization group and 
cluster approximation approaches agree for the most part. However, 
as mentioned, the agreement is not perfect, and is in any case limited 
to one particular model (the Farey fraction spin chain). Therefore 
further work on these models seems called for. 
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APPENDIX A: ASYMPTOTIC ANALYSIS OF THE CLUSTER 
PARTITION FUNCTION 

Here, we present the derivation of the asymptotic behavior (|38fl of the 
cluster partition function A(z,/3) near (z,[3) = (l,/? c ), where the critical 
temperature f3 c is determined by A(1,/3 C ) = 1. 

We recall that the cluster partition function is given by 

oo 

A(z,/3) = 5>"e-*». (Al) 

n=l 

where in our case 

e„ = ilog[»(n+l)], (A2) 
and A(l, 2) = 1 implies /3 C = 2. For (3 = (3 C one obtains 

A{z, &) = 1 + — log(l - z) . (A3) 

z 

For P > 1 and |z| < 1, A(z, /3) is an analytic function in j3 with coefficients 
depending on z. In particular, expanding around f3 — f3 c , we get 

A(z, /3) = A(z, C ) + V ([n(n + l)]<A-«/A> _ i) (A4 ) 

^-^ n (n + 1) V / 

n— 1 v ' 

= A(z, /3 C ) + g iog[„(„ + l)] + 0((/3 c - /3)2) 

p c n^n -t- j. j 

uniformly in \z\ < 1. In particular, letting z approach one from below, the 
coefficient of the linear term in (/3 C — (3) changes continuously with z and 
we arrive at 

A(z, (3) = A{z, f3 c ) + ^^(C + o((l - zf) + 0((/3 c - (3) 2 ) (A5) 

with C = J2n=i l °n(llt) )] = 2.046277452855878591 . . . (the sum can easily 
be evaluated numerically using the Euler-MacLaurin formula) . Introducing 
the reduced temperature t — 1 — (3/(3 c , this implies 

A(z,/3)~l + C< + (l-z)log(l-z) , (A6) 

which is our desired equation (|38|) . 

We conclude this appendix with a few generalizing remarks. For our 
purposes it was sufficient to work directly with l|22|) . but we would like to 
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point out that it is possible to perform a more thorough analysis. Observing 
that 

—_L-_= [°° K(P,s)e- ns ds, 
[n(n + l)f/ 2 J Q 

where 

ms) = ^ ) I i ,- i)/2 ( S /2)e-^-^, 

leads to 

00 7 n poo 

A(z,/3) = V = / #0M- da, 

an integral representation which is a different starting point for an asymp- 
totic analysis. In particular, one recognizes directly that the large-n asymp- 
totics of the cluster energies e n is related to the small- s expansion of the 
integral kernel K(/3, s), which in turn determines the singular behavior of 
A(z,/3). The argument can therefore also be extended to more general e n . 
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